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ABSTRACT 


Fractional and fractal derivatives are both generalizations of the usual derivatives 
that consider derivatives of non-integer orders. Interest in these generalizations 
has been triggered by a resurgence of clamor to develop a mathematical tool 
to describe "roughness" in the spirit of Mandelbrot's (1967) Fractal Geometry. 
Fractional derivatives take the analytic approach towards developing a rational 
order derivative while fractal derivatives follow a more concrete, albeit geometric 
approach to the same end. Since both approaches alleged to extend whole 
derivatives to rational derivatives, it is not surprising that confusion will arise 
over which generalization to use in practice. This paper attempts to highlight 
the connection between the various generalizations to fractional and fractal 
derivatives with the end-in-view of making these concepts useful in various 
Physics applications and to resolve some of the confusion that arise out of the 
fundamental philosophical differences in the derivation of fractional derivatives 
(non-local concept) and fractal derivatives (local concept). 

Keywords: fractional derivative, fractal derivative, fractional differential operators, fractal analysis 


I. NON-INTEGER ORDER OF DIFFERENTIATION 

L'Hospital (1695) wrote a letter to Leibnitz 
raising the possibility of a half derivative of the 
function f(x) = x to which the latter responded 
that "this could lead to a paradox for which 
the consequences may be useful in the future". 
Leibnitz realized ab initio that there are various 
ways of generalizing the discrete operation of 


differentiation to the continuous case which lead 
to non-equivalent results. To date, there are two 
(2) main streams of generalizations pursued 
by researchers: the purely analytical stream 
(Fractional Calculus) and the physical-geometric 
stream (Fractal Calculus) (West, Bologna & 
Grigolini, 2003; Miller & Ross, 1993; Hilfer, 2000). 
To illustrate the differences in results for the two 
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approaches, the half derivative of f(x] = x are: 

dlf (x) 2-fx 

1. dx*-/ 12 * -fir (fractional derivative] 

d/O) _ ,r- 

2. = — x (fractal derivative with 
scaled x] 

df 1 * 1 

3. diK i/a — 1 (fractal derivative with both 
variables scaled] 

The half fractional derivative of f(x] does 
not represent a local operation i.e. it does not 
represent "slope" at a point, but differentiation 
over a range of values of x. On the other hand, 
the half fractal derivative of the same function 
re-scales either the independent variable or 

■ y l-! 1 " 

both - . ) and applies the usual differentiation 
dx L ‘" 

on the re-scaled values. The latter approach is 
greatly influenced by the work of Mandelbrot 
(1967] in Fractal Geometry where description 
of the "roughness" of physical objects is the main 
focus. 

Recent developments in Fractional Calculus 
have demonstrated the intrinsic mathematical 
elegance of a purely analytic approach. For 
instance, the generalized differential operator 
defined by: 

4. 0'" [/(*)] = f ix) 

where B s f(x) — fix — e), yields the usual 
differentiation for n = 1 and the usual integration 
for n = -1. Such elegant results pave the way 
for a re-examination of the various concepts in 
classical Calculus viz. the Fundamental Theorem 
of Calculus. Meanwhile, developments in Fractal 
Calculus have demonstrated a strong inclination 
toward practical applications of the concept of a 
real order derivatives in physics and the physical 
sciences. Chen et al. (2003] used Fractal Calculus 
in investigating anomalous diffusion; Pant et al. 
(2012] considered Fractal Calculus in the study of 
diffusivity in porous media; Danijela et al. (2013] 


investigated Columbeau-generated stochastic 
processes through fractal derivatives. 

In a 2012 paper, Xiao-Jun Yang (2012] used 
the term Local Fractional Calculus (LFC] to 
denote Fractal Calculus for the development of a 
Local Fractional Fourier Analysis. He defines local 
fractional calculus as "dealing with everywhere 
continuous but nowhere differentiable functions 
in fractal spaces". We note that Yang's paper 
highlights the key difference between Fractional 
Calculus, on one hand, and Fractal Calculus on 
the other. The former is a non-local concept while 
the latter is a local one. In Fractional Calculus, we 
speak of differentiation over a range of values 
of x while in Fractal Calculus, we deal with 
differentiation at a point x. 

As a relatively young field of inquiry, there is 
bound to be a lot of confusion over the concepts 
of "fractional" and "fractal" derivatives because 
both deal with non-integer derivatives. The first 
concept is a purely analytic one while the latter has 
a strong geometric flavor. In fact, we propose to 
use the term "scaled calculus" for fractal calculus 
in order to reduce confusion with fractional 
calculus. The paper is organized as follows: 
Sections II and III deal with the development 
of fractional derivatives; Section IV derives 
fractal derivatives; Section V concludes with a 
series of simple applications of fractal calculus 
in Physics. With this paper, it is hoped that: (1] 
the differences, both in philosophy and use, of 
fractional and fractal derivatives are clarified, and 
(2] research interest in both developments will be 
further stimulated. 

II. GENERALIZATION TO FRACTIONAL 
DERIVATIVES BASED ON EXPONENTIAL AND 
POLYNOMIAL FUNCTIONS 

Consider the exponential function: 

1. f(x] = e ax 
whose kth derivative is: 














Maglasang, G., Ontolan, J. M., Borres, M. S. and Patac, A. V. 


221 


2 . 


dx* 


^-e ax = a k e ax 


A natural way to generalize to rational 
derivatives is to set k = q £ R + to yield: 


3. 


dx *? 


e- Q * = a? e a:i ' 


Any function that can be expressed as a linear 
sum of exponential functions can be differentiated 
in rational orders. In particular, since: 

, , 

4. cosU) = --- 

It follows that: 

5. ^ (COS<>)' = COS(> f ;?). 

The last result indicates that the fractional 
differential operator shifts the phase of the 
trigonometric function in fractional proportion. 
The same statement can be said of the sine function. 
The approach using exponential functions as a 
way to generalize fractional derivatives appears 
to be satisfactory yet, L'Hospital wanted the "half 
derivative" of f(_x) = x. There is no apparent 
representation of f(x) as an exponential function 
taking the form of (1). Note that f(x] = e log(x] is 
not the same as (1). The other way is to find the 
Fourier representation of the function (which is 
not defined over the entire interval] over some 
finite interval. However, ambiguities cannot be 
avoided here because of the choice of the finite 
interval. In other words, the fractional derivative 
obtained in this manner is neither unique nor 
local (i.e. whole derivatives are both unique and 
local]. 


Instead of beginning with an exponential 
function, we can attempt to build fractional 
derivatives by considering the polynomial 
function: 

6. fix) = x m 


whose kth derivative is : 

ftk 

7. = m(m — l)(m — 2) ,..(m — k + 

The Gamma function or generalized factorial 
can be used to make the result above applicable 
to any positive rational order of derivative. The 
Gamma function is : 


8. r(ff) = J c .^ LT 1 e v dx = {a — l)r(or — 1) ; 

which allows us to write (7] as: 

9 .4*- = .2=«L i -h. 

r(m+l-fe) 

with the recurrence relation T (n) = (n — l) I and 
the reflection relation r(x) P(l — x) — —-— 

siftTZX 


Replacing k with q £ Q+, we obtain the desired 
generalization. In response to Leibnitz, the half 
derivative of x now becomes: 



10 . 





1/2 


since P 




2 


While we are used to thinking of the 
derivative of a constant to be equal to zero, the 
half-derivative of a constant is not zero but: 



which is not defined at x = 0. 


We point out that combining the two 
approaches lead to some curious results. For 
instance, we know that: 





?a! 


but 






13. dar 1 .-' 3 


e x — 


gX — gX 




dx^ 


(e x ) - 


d L/a r n 

Li _ OG B **• 


dxy 
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Results such as (13] are not rare and for which 
reason, Leibnitz was right in pointing out that 
generalization to rational order derivatives 
indeed lead to paradoxes. 

III. GENERALIZATION TO FRACTIONAL 
DERIVATIVES BASED ON A LIMIT DEFINITION 

We next examine yet another approach that 
is more basic than the derivations of fractional 
derivatives based on the exponential and 
polynomial functions. The definition of a positive 
integer derivative is based on: 

14. f (fM) = 

Li. .r. £ 


For n = 1, we obtain the usual differential operator: 


17. D(fix ) ) = lim £ _»Q 


/UWO-s} 


while for n = -1, we make use of the geometric 
series representation: 


18. 


Ci-^} 


= e[1 + B s + ... -F 3 S + ...] 


This gives: 


19. zr 1 [/(*)]= Um^(^y^)/00 = 


E [ fix) -f fix — z) -F fix — 2z) -f ■■ ■ fix — kz) -f ... ] 


which can be repeated n times to yield: 

(/(*)) = Ey= 0 (-1) ; TlCj f(x -jz) 

where nC. is the number of combinations of n 

j 

objects taken j at a time. In order to generalize 
this formula when n is not a positive integer, we 
need to generalize in two directions: the binomial 
coefficients and the upper limit of the summation. 


Equation (16] gives the differential operator 
while (19] yields the integration operator in the 
limit as £ tends to zero. An algebraic expression 
for the Fundamental Theorem in Calculus can be 
obtained from (19] and (9] namely, that: 




dx k 



IT ( mil) 

P7i Ik = n -v ^ then for 
A U X , men ror 


k=-l, D 1 (x m ) is the integral of 


Consider a function f(x] shown below : 
Figure 1. Function f(x) 



Given £ > 0, we move five equi-steps to the left to 
obtain the points x - 5e, x - 4e,..., x - £.. In general, 
let k be the number of such points , then x = k£. 
Define the £- backward shift operator as: B £ (f(x] 
= f(x-e]. The general differential operator D is 
defined by: 

16. D"[/(*)] = f(x) 


For example, if ^("v) = th en 0 1 (:V 2 ) = 

-■ x 2-(-i) _ ri x 3 _ 1 v-3 

r(2+i-(~i)) 3! 3 ‘ 

The upper limit of the summation in Equation 
(15] becomes: 


Fr*] 


(/W) = X 


J=0 

where x c is usually zero, or 



We then proceed to define the binomial 
coefficients in a way that does not require integer 
values of n. Note that: 







which suggests that an alternative way of writing 
factorials for non-integer arguments be used. This 
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is facilitated by using the Gamma function [8]. 
Doing so, we obtain: 


X—X-r 


d q 


dx a - 


{fix)) = j (-1 V 


;-Q 


Hq + i) 
j' rfa + i-;) 


fix - /£■) 


The general derivative depends on the value 
of the function f over the range from x Q to x and 
is therefore a non-local operation. We can see 
this from the fact that in the factor f(x-j£] in the 
summation, the argument ranges from x to zero 
as j ranges from j = 0 to j =[(x — xO)/e]- 


Using (20], L'Hospital's query to Leibnitz yields: 


21 . 



A' 1 / 2 

V-ji 


which is identical with (8]. 


We therefore see that the generalized derivative 
obtained by the polynomial approach and 
that which is obtained by the limit approach 
give identical results. The utility of using the 
generalized limit approach over the polynomial 
approach is that the generalized limit approach 
allows us to find rational order derivatives even if 
the function is not expressible as a power series. 
Moreover, the polynomial approach does not give 
us any indication about the non-locality of the 
generalized derivative, that is, the dependence 
on the function over a range rather than on just 
a single point. 


IV. RECONCILING THE EXPONENTIAL 
APPROACH AND THE GENERALIZED LIMIT 
APPROACH TO FRACTIONAL DERIVATIVES 

In Section II, we pointed out the conflicting 
results when we took the half-derivative of 
e x using the exponential approach and the 
polynomial approach (or the generalized limit 
approach]. The derivative of this function 
(including half derivatives] should yield the same 
function. Moreover, the function approaches 1 as 
x tends to zero. Yet, the half derivative of e x is: 


x—x, 


1/2 ^ V—i 

^ IT {f{x))= lim^ 2. (-1V jy 


T{3/2) 


}=a 


j'. r(3/2 -j) 


fix - JB ) 


1/--II-T , 2i ® 3 „ Ik 4, 

-==. 1 + 2x -I- x. --I- x 1 T »' 

-ff TTX - 3 15 105 


which goes to infinity as x goes to zero. How can 
we reconcile the two results? The answer lies 
on the fact that half-derivatives depend on the 
chosen range of differentiation. For instance, 

e x 15 

and the only way for the derivative of the right - 
hand side to equal to the argument of the definite 
integral on the left-hand side is when _*: c = — oo 
which we have assumed tacitly. In other words, 
ranges of differentiation are important when we 
talk of fractional derivatives (derivatives over a 
range of values rather than derivative at a single 
point]. If x c = —oo, then we find that the upper 
limit of (22] simply becomes oo and hence, the 
half-derivative in (21] becomes e x also. 



V. FRACTAL DIMENSIONS, ROUGHNESS, 
FRACTAL CALCULUS AND FRACTIONAL 
CALCULUS 

The motivation for developing Fractional 
Calculus stemmed from a purely analytic question 
of L'Hospital (1695]. Recent interest in fractal 
geometry and analysis, however, spawned yet 
another direction for generalizing the usual 
derivative. This direction is towards defining 
Fractal Derivatives as basis for describing data 
and geometric roughness. As such, we expect 
a relationship between fractal derivatives and 
fractal dimensions. What is not clear yet is the 
relationship between fractal derivatives and 
fractional derivatives. 


Fractal Derivative 

Let f(t] be a function. We define the q th fractal 
derivative of f(t] with respect to a fractal measure 
t as: 


(23) 


/(0 

dr? 


= lim 





,q > 0 
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Note at once that while for fractional 
derivatives £ = 1 £ — t ± | for fractal derivatives 
we have — \^ q — t q |- This implies that 


(24) 




Comparing [24] and (20], we see that (24] 
cannot be expressed in terms of simple backshift 
operators. In particular, the half - derivative of 
m = £ becomes: 


(25] 



i 

which differs from (21] by a factor of A more 
general representation of a fractal derivative with 
respect to both f(t] and t is: 


e ,fl_ e q 


Hence for f(t] = t with r — q = 7 : 


(27] 


d/ 


l/E 


d tl 


- lim ^=—= 1 


which coincides with — = 1 , 

dt 


Observations: 


where T(.] is a functional of G (■ > ■ ) - 

We can compare the fractal derivative G(0 r 
of fit) =6" with its corresponding fractional 

derivative -4 For reference, we reproduce 
( 22 ] below: 

A (— i [e r ] = -4 (H + 2t5 + -tl + — £5 + — tl ... + 

\2J Vnr 4 lb 1GS 

-> Fractional Derivative 

The corresponding fractal derivative is: 


G{Q^)[e f . = 2U + 2ta -h f ....4 —>Fractal Ztenvatiw 

For £ > 0, the fractal derivative is generally 
greater than the fractional derivative. As t > 0 
, the fractional derivative tends to infinity 
while the fractal derivative tends to zero; both 
behaviors being inconsistent with the fractional 

f 

derivative of e ' using the exponential approach 
viz. f f 'v {e c ") = 1 as t -> 0, The 

behaviors of the derivatives A £), G(0, i) and 

ft (t) of e : are shown below 


df 

a] The fractal derivative d q is always a function 

d*f 

of the fractional derivative d{ _ q . 

b] To avoid confusion, we change notation for 


J JfT 


fractal derivatives. Instead of writing —, we 

itfl 

write G[r } q) for Equation (26]. Likewise, 

instead of writing d f° r the q th fractional 

derivative, we write A(q]. The ff (., . ) notation 
reflects the geometric nature of (26] while 
the A (.) notation highlights the analytic 
feature of the fractional derivative Equation 
( 20 ] 

c] Statement (a] can be written as: 

A(?) = r(G(Q, q)), 


Figure 2. e l 



















G(0,1/2) 
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Figure 3. half fractal derivative 


Figure 4. half fractional derivative 
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0 < t < .51 

1 

1 

2 

t 

1 

3 

0 

1 < t < 1.036 

i 

2 

t 

3 

t 

exp(t) 

A(l/2) 

G(0,l/2) 

t 

exp(t) 

A(l/2) 

G(0,l/2) 

0.00169 

1.0017 

13.77046 

0.082358 

1.00019 

2.7188 

2.502952 

4.66774 

0.00622 

1.0062 

7.242892 

0.158716 

1.00045 

2.7195 

2.503464 

4.6693 

0.00636 

1.0064 

7.164721 

0.160515 

1.00645 

2.7359 

2.515332 

4.705406 

0.02289 

1.0232 

3.901294 

0.309542 

1.01012 

2.7459 

2.522633 

4.727583 

0.0237 

1.024 

3.84009 

0.315223 

1.01249 

2.7524 

2.527365 

4.741942 

0.02547 

1.0258 

3.717008 

0.327352 

1.01273 

2.7531 

2.527845 

4.743398 

0.02599 

1.0263 

3.683348 

0.330845 

1.01772 

2.7669 

2.537853 

4.773732 

0.02788 

1.0283 

3.569347 

0.343301 

1.02446 

2.7856 

2.551465 

4.814913 

0.03151 

1.032 

3.381066 

0.366268 

1.02761 

2.7944 

2.557863 

4.834241 

0.50395 

1.6552 

1.80922 

2.225664 

1.03589 

2.8176 

2.574793 

4.885299 

1.5 < t <1.55 




2 < t < 2.05 




T 

exp(t) 

A(l/2) 

G(0,l/2) 

t 

exp(t) 

A(l/2) 

G(0,l/2) 

1.50904 

4.5224 

3.83446 

8.503707 

2.00068 

7.3941 

5.86979 

14.15103 

1.5178 

4.5622 

3.863687 

8.585542 

2.01024 

7.4651 

5.918051 

14.28454 

1.51835 

4.5647 

3.86553 

8.590701 

2.01281 

7.4843 

5.931086 

14.3206 

1.51986 

4.5716 

3.870594 

8.604875 

2.0134 

7.4887 

5.934082 

14.32889 

1.52179 

4.5804 

3.877077 

8.623018 

2.01531 

7.503 

5.943791 

14.35575 

1.52312 

4.5865 

3.881551 

8.635538 

2.03163 

7.6265 

6.02734 

14.58689 

1.52507 

4.5955 

3.88812 

8.653921 

2.03406 

7.6451 

6.03987 

14.62156 

1.52573 

4.5985 

3.890347 

8.660149 

2.03733 

7.6701 

6.05677 

14.66831 

1.52706 

4.6046 

3.894837 

8.672711 

2.04045 

7.6941 

6.072934 

14.71304 

1.53064 

4.6211 

3.90695 

8.706595 

2.04153 

7.7024 

6.078538 

14.72854 
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In the neighborhood of zero, the analytic 
fractional derivative is larger than the geometric 
fractal derivative. As we move away from zero 
(somewhere near t = 0.50], the geometric fractal 
derivative dominates both the analytic fractional 
derivative and the fractional derivative obtained 
through the exponential approach. There appears 
to be some simple relationship that governs the 
behavior of the exponential fractional derivative 
(et) and the geometric fractal derivative. We 
performed a regression analysis using y = et and x 
= G(0,l/2]. The results are shown below: 


Figure 5. Graph of x = G(0,l/2) versus y = exp(t) 



The regression equation is 
exp(t) = 0.198 + 0.539 G(0,l/2) 

Predictor Coef SE Coef T P 

Constant 0.19798 0.01841 10.75 0.000 

G(0,l/2] 0.539015 0.001241 434.33 0.000 

S = 0.3696 R-Sq = 99.5% R-Sq(adj) = 99.5% 


That is, the exponential fractional derivative is 
about one-half of the geometric fractal derivative. 
The two plots of the other possible relationships 
between the three half-derivatives are shown 
below: 

Figure 6. Plot of x = A(l/2) and y = exp(t) 


Figure 7. Plot of x = A(l/2) and y = G(0,l/2) 



i-1-1-1-1-1-1-r 

0 2 £■ 6 3 *0 *2 "4 

A(1'2) 


Note that Figures 5 and 6 are almost identical. This 
observation simply comes from the fact that the 
geometric fractal derivative (G(0,l/2]] and the 
fractional derivative obtained by the exponential 
approach (e 1 ] are linearly related. 


VI. SOME APPLICATION OF FRACTAL 
DERIVATIVES 

Physical laws of diffusion that are based on 
a Euclidean medium cannot apply if in fact the 
media exhibit fractal properties. Some physical 
realizations of fractal media include water 
turbulence porous media and aquifers. In the 
latter case, the physical concepts of space and 
time have to be rescaled to, ( x r , t q ] where r, q 
are positive real numbers, in order to conform 
to the fractal nature of space-time. In this fractal 
space-time system, the notion of a velocity can be 
redefined as: 

28. G[r, q)S(x,t) - ^ 

because the traditional definition of velocity is 
meaningless in a non-differentiable fractal space- 
time. 



Fractal Linear Motion 

Linear motions involving the concepts of 
velocity and acceleration are fundamental in 
Physics. The usual distance traversed by a free- 
falling body after time t is: 

29. x(tj = -gt 2 
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where e is assumed constant. Suppose that the 
space *0) is fractal of dimension r , then the 
velocity of the body at any constant t, becomes 


30. 




£± t 


= lim 


r 

1- yl 




t-ti 


— , 2 r ~ 1 
- 2 r-±g ti 

Note that if the fractal dimension of Jt(t) is 

r = ljWe obtain the usual velocity' as: 


equation obtained by transformation of variables. 


Fourier's Law of Heat Conduction 

As a third example, consider Fourier's 
law (1822] (see Fourier, 1955] which states 
that "thermal conduction results into heat flux 
which is proportional to the magnitude of the 
temperature gradient and opposite to it in sign". 
Thus, for a one-direction conduction process, this 
can be translated to: 



dT 



v = -gt 

However, if x(t) is fragmented with, say, 
r - 0.63;the velocity is: 

w M3 =-(2 a37 a:Ct).63)fl 0 - H t 1 a36 ) 

which is slower than gt, as one would expect 
if one were travelling along a rough media e.g 
variable air resistance. 


where c is the thermal conductivity constant, h x 
is the heat flux in W/m 2 in the positive x direction 
and dT/dx is the (negative] temperature gradient 
(K/m] in the direction of the heat flow (usually, 
from hot to cold regions]. Suppose now that the 
space (x] is fractal of dimension r while T is of 
fractal dimension q, then Fourier's law translates 
to: 

. dT A 

32. h V = —c —- 


Fick's Second Law Anomalous Motion 

Consider the anomalous transport-diffusion 
process which is an alternative to Fick's Second 
Law given by: 


d-u(j£7t) 




d fdu(x f t) 


dx' 


—co < x < co ; u (x f 0.) = 8(x) 


where 5(x] is the Dirac delta function, and 0 < q < 
2, 0 < r < 1. The fundamental solution to (32] can 
be obtained by the transformation t* = t q and x* = 
x r which transform (32] into the usual diffusion 
equation. The solution, of course, is given by: 

33 -“ fet)= 

which is the stretched normal (Gaussian] 
distribution. From this example, we learn that 
in order to deal with fractal space-time, it is 
necessary to know the fractal dimensions r and 
q of the space (x] and the time (t] respectively. 
The corresponding fractal derivatives can be 
calculated and the solution to a fractal differential 


whose solution is obtained by the transformation 
x* = x r and T* = T q . Equation (35] is the anomalous 
heat conduction equation. 

Other examples may be constructed where 
the use of fractal derivatives is examplified. The 
reader may wish to refer to Darcy's law . Darcy's 
law relates the instantaneous discharge rate 
through a porous medium and the viscosity of the 
fluid with pressure drop over a given distance. 

VII CONCLUSION 

While we have found a linear relationship 
between fractal derivatives (G(0,q]] and fractional 
derivatives through the exponential definition, 
much is still to be done in terms of relating 
analytic fractional derivatives (A(q]] and fractal 
derivatives (G(0,q]]. A relationship between A(q] 
and G(0,q], if found, would be useful in terms of 
developing an entire theory of fractal calculus (as 
opposed to fractional calculus]. A fractal calculus 
theory would, in turn, be useful in re-analyzing 
most of the laws in Physics e.g. quantum theory. 
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The key to such an examination lies on the 
premise of a local versus non-local derivatives. 
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